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We show that the observational data recently provided by 
Giovanelli et al. (1996 a, b) and discussed by BahcaU and 
Oh (1996) concerning the velocity distribution of clusters of 
galaxies can be naturally fitted by a statistical distribution 
which generalizes the Maxwell-Boltzmann one (herein recov- 
ered for the entropic index q = 1). Indeed a recent general- 
ization of the Boltzmann-Gibbs thermostatistics suggests for 
this problem that the probability function is given, within a 
simple phenomenological model, by 



P{>v) = 



dv[l~{l-q){v/vof 
dv[l-[l-q)(v/v,Y 
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g)^^/^ ifg<l 
if g > 1. 



A remarkably good fitting with the data is obtained for q 



0.23 



+0.07 



and «o ~ 490 ± 5 km 



The motions of clusters and group of galaxies have 
been studied by BahcaU et al. (1994). The sample of 
clusters of galaxies velocities recently observed by Gio- 
vanelli et al. (1996 a, b) based on TuUy-Fisher distances 
of Sc galaxies has been examined and compared with 
model expectations by BahcaU and Oh (1996). In con- 
trast with previous analysis of the data, the actual ob- 
served velocity function does not show a tail of high ve- 
locity clusters (clusters with velocities greater than ~ 600 
km s^^ can be found with a probability less than five per 
cent). BahcaU and Oh (1996) determine the cluster ve- 
locity function using simulations based on four different 
cosmological models: 

(a) CDM (cold dark matter), = dimensionless matter 
density = 0.3; 

(b) FBI (primeval baryonic isocurvature) , fl — 0.3; 

(c) CDM, Q = 1.0; 

(d) HDM (hot dark matter), Q = 1.0. 

These authors conclude that model (a) and marginally 
model (b), in spite of the high energy tail, are consis- 
tent with the observed data (see Fig. 1), while models 
(c) and (d) show a too large high velocity tail to be in 
agreement with data. However, if we look in detail at the 
observations and simulations curves, we can see that not 
even the CDM, il — 0.3 model is completely satisfactory, 
particularly in the region of velocities larger than 350 km 
s~^. In fact, none of the four theoretical curves is fully 



satisfactory. 

We want to show in this letter that to understand 
the Giovanelli et al. (1996 a, b) observational data is 
not simply a problem of which cosmological model is 
used, but it also depends on which statistical mechan- 
ics is used. In fact, if we adopt the recently introduced 
non-extensive statistics (Tsallis 1988, Curado and Tsal- 
lis 1991) to calculate the probability distribution function 
of cluster peculiar velocities, the agreement with exper- 
imental results is very satisfactory, particularly in fore- 
seing that the number of galaxy clusters with velocities 
greater than 600 km s""'^ is practically zero. The non- 
extensive statistics, among several other applications, has 
been recently applied in calculating matter distribution 
of self-gravitating systems (Plastino and Flastino 1993, 
Aly 1993, Boghosian 1996), turbulence (Boghosian 1996), 
Levy- and correlated-like anomalous diffusions (Zanette 
and Alemany 1995, Tsallis et al. 1995, Tsallis and Buk- 
man 1996), solar neutrino fluxes (Kaniadakis et al. 1996), 
a cosmological model (Hamity and Barraco 1996), as well 
as in developing a dynamic linear response theory (Ra- 
jagopal 1996). 

We firstly give a brief account of the non-extensive 
statistics. Then we calculate by means of it the probabil- 
ity of finding galaxy clusters of certain velocity v. Results 
are shown in the two Figures and discussed; conclusions 
are outlined. 

In order to cover non-extensive systems (long-range mi- 
croscopy memory, long range forces, fractal space time) 
the following generalized entropy has been proposed: 
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where fc is a positive constant. Optimization of Sq yields, 
for the canonical ensemble. 
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and, when q ^ 1, the Boltzmann-Gibbs result {pi cx 
e^^'^*) is recovered. Let us stress that, for g < \, pi 
vanishes for [1 — (1 — q)(3ei\ < 0. Within this formalism 
the observational data are to be identified with {A)q — 
pfAi where {Ai} are the eigenvalues of an arbitrary 
observable A. 
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The apphcation of this thermostatistics to stellar sys- 
tems has shown that sensible distribution functions (stel- 
lar polytropes) can be derived if 9 G (—00, 7/9). Distribu- 
tions with q < 1 give rise to the spatial cutoff of the mass 
distribution and the finite mass of the stellar polytrope 
(Plastino and Plastino 1993, Aly 1993 and Boghosian 
1996). Observed distributions of pure-electron plasma 
with two dimensional turbulence are well described with 
q = 0.5 (Boghosian 1996). The logistic map at its thresh- 
old to chaos can be described using q = 0.24 (Tsallis et 
at 1996). Solar neutrino fluxes can be understood using 
non-extensive statistics for the central core plasma and 
values of q slightly below one: q = 0.997 (ion plasma ) 
and q = 0.976 (electron plasma) (Kaniadakis et al. 1996, 
Quarati et al. 1996). 

We calculate now the probability P{> v) of finding 
one-dimensional peculiar velocity function of clusters of 
galaxies (relative to a comoving cosmic frame) greater 
than V. It is given by 



P{> v) 
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where vq is a phenomenological characteristic velocity 
which incorporates all the averages over cluster config- 
urations taking into account their gravitational interac- 
tions. 

We find the best fit of the observed data and report it 
on Fig.l. The entropic parameter is g = 0.23^0 05 and 
Vq ~ 490 ± 5 km s~^. This value of q, very different 
from one, reflects the fact that the long-range gravita- 
tional forces are, of course, in the case of galaxy clus- 
ters, very important. Moreover, by means of the relation 
n = polytropic index = 3/2 + q/{l — q) (Boghosian 
1996), we may derive that the system is characterized by 
n « 7/4. The fact that this value lies in the range of 
the typically observed (and theoretically allowed: Plas- 
tino and Plastino 1993, Aly 1993, Boghosian 1996) values 
constitutes a further, and independent, confirmation of 
the possible relevance of the present proposal. 
Peculiar of the non-extensive statistics is that our distri- 
bution is cut at V > 520 km/s, in perfect consistence with 
the observed data. In Fig. 2 we show P(> v) of Eq.(|) 
for six different values of q: q = 1 (implicitely adopted 
by Bahcall and Oh 1996) is the Maxwell-Boltzmann dis- 
tribution, (7 = is a straight line, the q < curves show 
clearly a pronounced cut tail. 

Our results indicate that in order to distinguish which, 
among the four models judiciously used by Bahcall and 
Oh (1996) (each of them using several fitting parame- 
ters) , is the most appropriate to describe the velocity dis- 
tribution, it would be necessary to repeat the simulations 
within the non-extensive statistics with q = 0.23lg'Q5. 



The best value of vq can be selected and, possibly, one 
could obtain information on the quantity fi. Finally, we 
note that the efficiency coming from the modification of 
the statistics is greater and more important than the effi- 
ciency obtained from the modification of the model (pre- 
cisely the same occured for experiments done in pure- 
electron plasma turbulence: see Huang and DriscoU 1994 
and Boghosian 1996). In fact, our fit of the experimental 
results, based on the non-extensive statistics, has been 
derived within a classical ideal gas picture. Incidentally, 
we notice that the value q = 0.23^Q g5 is very close to the 
value q = 0.24, which is appropriate to chaotic dynamics 
(Tsallis et al. 1996). It would be interesting to investi- 
gate if this approximate equality between the two above 
figures is just a numerical coincidence or is due to some 
fundamental reasons. 
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Figure Captions 

Fig. 1. - The observed data of the clusters velocity 
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function arc reported together with the four curves cor- 
responding to the four different models of P{> v) elab- 
orated by Bahcall and Oh (see their Fig. lb). Our fit 
of the observed CVF, based on non-extensive statistics 
{q = 0.23 and vq = 490 km s^^), is also shown. 
Fig. 2. - The function P(> x) versus x = v/vq is re- 
ported for six different values of the cntropic parameter 
q (—5, —1, 0, 0.5, 1, 2). The case g = 1 is the well known 
Maxwell-Boltzmann distribution. We note that the curve 
5 = is a straight line. Curves with negative values of q 
have an abruptly cutted tail. 
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